Resonances in a dilute gas of magnons and metamagnetism of isotropic frustrated 

ferromagnetic spin chains 
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We show that spin-S chains with SU(2)-symmetric, ferromagnetic nearest-neighbor and frustrating 
antiferromagnetic next-nearest-neighbor exchange interactions exhibit metamagnetic behavior under 
the influence of an external magnetic field for small S, in the form of a first-order transition to the 
fully polarized state. The corresponding magnetization jump increases gradually starting from 
an 5-dependent critical value of exchange couplings and takes a maximum in the vicinity of a 
ferromagnetic Lifshitz point. The metamagnetism results from resonances in the dilute magnon gas 
caused by an interplay between quantum fluctuations and frustration. 



I. INTRODUCTION 

Quantum spin systems at low temperatures share 
many of the macroscopic quantum behaviors with sys- 
tems of bosons such as Bosc-Einstein condensation, su- 
perfluidity (see Ref. [l] and references therein) , or macro- 
scopic quantum tunneling^ and may indeed be viewed 
as quantum simulators of interacting bosons X Spin sys- 
tems with frustration, in particular, realize exotic phases 
of strongly correlated bosons, such as spin liquids 3 - or 
supcrsolids.— Experimentally, such systems can be stud- 
ied both in low-dimensional quantum magnets (see, e.g., 
Ref. H) and in ultra-cold atomic gases^ In the latter 
case, interest in the many-body physics has been excited 
by the extraordinary control over interactions between 
bosons via Feshbach resonances^ These can, in particu- 
lar, be used to tune interactions from repulsive to attrac- 
tive. In the attractive regime, the collapse of an ultra- 
cold gas of bosons was observed in experiments^ In this 
work, we discuss a mechanism by which the same can be 
achieved in spin systems. 

In the large-S* limit, spins map onto bosons with a 
finite but large Hilbert space, justifying a description in 
terms of soft-core bosons. Yet, in many cases, similarities 
between spins and soft-core bosons may even exist for 
5 = 1/2. In three dimensions (3D), systems of spins 
and bosons resemble each other the more the smaller the 
density of bosons is, whereas ID spin-S 1 antiferromagncts 
close to saturation behave as spinlcss fermions, or hard- 
core bosons. In spin systems, the external magnetic field 
h tunes the density of magnons. The limit of a dilute gas 
of magnons is then realized as the fully polarized state 
(the vacuum of magnons) is approached from below. 

In this work, we argue that resonances can play a cru- 
cial role in frustrated quantum spin systems largely deter- 
mining their low-energy behavior in an external magnetic 
field. We consider the frustrated ferromagnetic (FM) 
spin- S chain and show that upon changing system pa- 
rameters such as the coupling constants or h, one can 



tune the effective interaction between magnons from re- 
pulsive to attractive by exploiting the existence of res- 
onances. As a main result, we demonstrate that, close 
to resonances, where the scattering length is much larger 
than the lattice constant, and in the case of attractive 
effective interactions, the intrinsic hard-coreness of spins 
does not play a significant role in the limit of a dilute gas 
of magnons. A behavior resembling the collapse of attrac- 
tively interacting boson s 9 ' 10 therefore exists in such spin 
systems close to their fully polarized state. The thermo- 
dynamic instability of collapsed states causes jumps in 
the magnetization curve just below saturation. This has 
to be contrasted with the scattering length being of the 
order of a few lattice sites. In that case, which is real- 
ized for spin 1/2, the formation of mutually repulsive, 
multi-magnon bound states can be observed (see, e.g., 
Refs. [lll-lll). 

Our work demonstrates that the mapping of a purely 
ID spin system close to saturation to an effective theory 
of a dilute Bose gas properly accounts for the physics of 
the model. This is accomplished by connecting the scat- 
tering vertices of the microscopic lattice model with the 
coupling constants of the effective theory. Furthermore, 
despite the purely ID nature of our problem, a 1/S ex- 
pansion is a valuable tool and yields the correct physics. 

Concretely, we study the following system: 
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J Si • Si. 



J Si ■ Si. 



hSf 



(1) 



Si = (Sf,Sf,S?) is a spin-S* operator acting on site i 
and L is the number of sites. J < is the FM, nearest- 
neighbor exchange interaction and J' = 1 is the antifer- 
romagnetic (AFM), next-nearest-neighbor exchange in- 
teraction setting the energy scale. In the absence of an 
external field h, Hq has a FM ground state for J < —4 
for all S (see Ref. Il6f) . J = —4 is a ferromagnetic Lifshitz 
point where the quadratic term in the dispersion of the 
magnons vanishes. 
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Systems with competing FM and AFM interactions are 
of timely interest ji 7 - in particular, the spin- 1/2 version of 
Eq. ([1]),— ~ 14 i 18 ~ — motivated by the experimental real- 
izations in, e.g., LiCuV0 4 (Refs. [HUH) and Li 2 ZrCu0 4 
(Ref. HI). The ID case of J < and 5 > 1/2 is largely 
unexplored (for J > and 5 > 1/2, see Refs .[l9l and Im). 

We are mainly interested in the region — 4 < J < 
and magnetization M = S Z /{SL) (S z = J2t( S t)) closc 
to saturation M = 1, for general spin 5. Wc will proceed 
in three steps: First, in Sec. [TTl we discuss the solution 
of the two-magnon problem and introduce the scatter- 
ing length. Second, in Sec. IIII( we map the low-energy 
limit of Eq. ([I]), close to saturation, to a dilute ID gas of 
two species of bosons interacting via an effective short- 
range interactioni 25 i 26 We then calculate the interaction 
vertices in this effective theory using a 1/5* expansion. Fi- 
nally, in Sec. IIV1 we compare the analytical results with 
exact numerical ones using the density matrix rcnormal- 
ization group (DMRG) metho d 27 ' 28 and exact diagonal- 
ization (ED). We put a particular focus on the case of 
5 = 1. A summary of our results is presented in Sec- 
tion [V] while technical details of the mapping to a di- 
lute gas and of the 1 / S expansion are given in Appendix 
|A"1 A comparison of DMRG results for open boundary 
conditions vs. results for periodic boundary conditions is 
shown in Appendix [Bj 



II. TWO-MAGNON PROBLEM 



A. Solution of the two-magnon problem 



This leads to the recurrence relations 

5 



£IqCq 



y/S(2S-l) 



(25 - I) 3 / 2 
(fio-J)Ci = 1 S3/ J (iC + (iC 2 

+(2(0, + C 3 ) 

(2S - 1) 3/2 

(no-l)Ca = S 3/2 C 2 Q) + C2C4 

+Ci(^i + c 3 ) 

OoCV = Cl (C r +i + C r _i) 

+C2 {C r+ 2 + C r - 2 ) , forr>3, (5) 

where Cl = 25Jcos (K/2), (2 = 2Scos{K). When 
|*2Af) is a bound state, O = E b -4S(1 + J 2 /8) where E b 
is the (negative) binding energy (defined as the bound- 
state energy minus the energy of the minimum of the 
two-magnon scattering states) . 

The (unnormalizcd) two-magnon bound states for a 
given K are constructed with the ansatz 



(r>l), 



(6) 



which, inserted in Eq. ([5]), leads to a characteristic quar- 
tic equation for r > 3 



n z 2 -^{z 3 + z)-c 2 {z A + i) = {), 



(7) 



z being any of e _K± with Rc[/c±] > 0. The remaining 
unknown quantities Co, v and E b are determined from 
the remaining relations listed in Eq. ([5]). 



We now solve the interacting two-magnon problem, 
starting with the thermodynamic limit (for 5 = 1/2, see 
Ref. [290 . On a chain of finite length L with periodic 
boundary conditions the total momentum K is a good 
quantum number due to the translational invariance of 
the Hamiltonian Hs in Eq. (JTJ) . Thus, it is convenient to 
use a basis separating momentum subspaces 



\K,r) = ^e iK ^S r S r+r \F), 
1=1 



(2) 



where \F) is the fully polarized state, K = 2qn/L 
(q = 0, 1, ■■ • ,L— 1) and r is the relative distance of 
two magnons. The allowed values of r depend on S and 
the parity of L and q. For instance, in the case of S > 1/2 
and L even, r = 0, 1, .., L/2 — 1, (L/2) for q odd (even). 

We expand a general two-magnon state with momen- 
tum K into the (unnormalized) basis of Eq. @ as 



B. Scattering length in the lattice problem 

For S > 1/2, bound states with energies below the 
minimum of the two-magnon scattering continuum ex- 
ist only for K ~ ±2k c i and only in a finite window of 
couplings 



4 < J < J cr (S), 



(8) 



with S'-dependent critical values J cr (5 l ), as illustrated in 
Table |H The critical value J cr (S), which is J cr w —2.11 
for S — 1, quickly approaches J cr (S) ~ —4 with increas- 
ing S (see Tab. fl}. In fact, the 1/5 analysis to be pre- 
sented in Sec. IIIII suggests the existence of an S cr beyond 
which this window disappears completely. 

We define the scattering length of bound states, in the 
thermodynamic limit L — > 00, from their spatial extent 



(3) 



s 
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3/2 


2 


5/2 




2.11 (2.95) 


3.31 (3.42) 


3.68 (3.66) 


3.84 (3.80) 



and determine C r analytically by solving the two-magnon 
Schrodingcr equation 



Hs\^2M) — E2MV&2M) 



(4) 



TABLE I: Critical exchange couplings J cr (S) for the existence 
of metamagnetism in Eq. (TTJ derived from solving the two- 
magnon problem (values in parenthesis: Results from the 1/5 
expansion of Sec. IIII|) . 
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states are not formed below the minimum of the scat- 
tering continuum). In the attractive regime as > 0, the 
scattering lengths obtained from both approaches [i.e., 
Eq. (JHJ) and Eq. (p~5]) ] are in excellent agreement with 
each other. As a side note on terminology, we call at- 
tractive (repulsive) regime the one in which the effective 
interaction between magnons is attractive (repulsive). 

We can now generalize the procedure^ of mapping the 
antifcrromagnctic (unfrustrated) spin-S* chain close to 
saturation onto the low-density limit of the Lieb-Liniger 
model with a coupling constant 



FIG. 1: (Color online) ID scattering length as /a (a: lattice 
spacing) for S = 1/2, 1,3/2 (solid, dashed, dot-dashed line). 



(in analogy to the continuum problem of particles inter- 
acting via a short-range, attractive potential): 



9o 



as 



mm{Re[K±]} 



(9) 



The binding energy takes its lowest value (i.e., the largest 
absolute value) for K = K* ~ ±2k c i = ±2 arccos (—J/4) 
and this quantity, with extremely high accuracy, is re- 
lated to the scattering length by 



E b {K*) 



where m is the one-magnon mass, 



5(4- J) (4+ J)' 



(10) 



(11) 



The relation Eq. (|10l) between the binding energy and 
the scattering length that holds for our microscopic lat- 
tice model is typical for a ID Bose gas in the continuum 
interacting via an attractive contact potential, the Lieb- 
Liniger model.— 

The scattering length as can also be determined from 
the scattering problem of two magnons for general S in 
the thermodynamic limit. We have solved this problem, 
with the momenta of the two magnons (which partici- 
pate in scattering) taken in the vicinity of the same dis- 
persion minimum, k\ = k c i + k and ki = k c i — k. From 
the asymptotic form of the two magnon scattering state 
wavefunction we extract the scattering phase shift Ss{k) 
for any S, 



lim C r ~ cos (rk + 6s(k)). 



(12) 



To extract the scattering length from the scattering phase 
shift we use the same relation as in a ID continuum model 
of particles interacting via a short-range potential 



as = lim 

fe->o 



cot (S s (k)) 



(13) 



This allows us to calculate the scattering length in the re- 
pulsive regime as < as well (when two-magnon bound 



mas 



(14) 



However, in our model the single-magnon dispersion 
has two minima. The effective theory thus will be a two- 
component (two species) Lieb-Liniger model. There are 
two types of low-energy scattering processes, hrst when 
momenta of two magnons are in the vicinity of the same 
dispersion minimum that we have presented above (in- 
traspecies scattering), and second, when the momenta 
k\ and ki of two magnons are in the vicinity of different 
minima of dispersion, i.e., k\ = k c i + k and ki = — k c i — k 
(interspecies scattering). For the latter case we can re- 
peat all steps presented above and extract another cou- 
pling constant go from the interspecies scattering length, 
as in analogy with Eq. (fl"4]l . We obtain that go > (im- 
plying that bound states with total momentum K = 
are never formed below the scattering continuum), and 
go > go for any S in the region —4 < J < 0. Since 
the relation g > go always holds, the relevant scatter- 
ing length at low energies is the intraspecies scattering 
length as- 

The scattering length, shown in Fig. [T] can be well 
described, for small S > 1/2, by a sum of two terms 
(resonances): as — A^/(4+ J) + Xg/(J cr (S) — J) [where 
Xg are numerical prefactors]. We emphasize that, for 
S > 1, the scattering length is in general much larger 
than the lattice spacing, as is evident from Fig. [1] For 
S = 1, as takes a minimum at J ~ —3.3 with as — 80a. 
Additionally, the emergence of bound states manifests 
itself by a diverging scattering length at J cr (S), where as 
changes its sign jumping from — oo to +oo. Thus, bound 
states are typically shallow, with a binding energy given 
by Eq. (|10p . In addition, the minima in their dispersion 
occur at incommensurate momenta K* . 

For S = 1/2, any J < induces a two-magnon bound 
state with total momentum K = 7T ) 18 ' 20 i.e., as > 
and there is no resonance at —4 < J < for S = 1/2. 
Hence, S — 1/2 is very different from the S > 1/2 case 
where bound states with K = tt are never below the two- 
magnon scattering continuum, and, as discussed above, 
a resonance exists for —4 < J < and 1 < S < S cr . 

In order to analyze the finite-size effects with respect 
to results in the thermodynamic limit, we have numer- 
ically diagonalized Hs in the basis given in Eq. © 
(see Ref . Qjl for details of the procedure) . We obtained 
the full spectrum, i.e., the scattering continuum and 
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oo 


-J cr {S=l) 


2.25 


2.16 


2.13 


2.11 



TABLE II: Finite-size dependence of critical values J cr (S) 
for the emergence of bound states below the minimum of the 
two-magnon continuum of scattering states, for 5 = 1. 



bound/antibound states if present, for selected values 
of J G [—4, 0] in systems with up to L = 4000 sites 
and several values of S. Table [TT| shows the numerical 
determination of J cr (S) for S — 1 and different sys- 
tem sizes. Although finite-size effects are apparent, a 
quadratic fit in 1/L to numerical data for J cr (S) extrap- 
olates to J cr ~ —2.11 in agreement with the result de- 
termined directly in the thermodynamic limit (see the 
preceding discussion and Table [I]). 



III. MAPPING OF THE SPIN HAMILTONIAN 
TO A DILUTE GAS OF BOSONS 



In this section, we describe our effective theory in the 
thermodynamic limit, for the case of a finite (though 
vanishingly small) density of magnons. The mapping 
to a dilute gas of bosons is motivated by the following 
observation: For S > 1/2, we have shown that ag is 
large. Hence in the dilute limit, we can safely neglect 
the hard-core constraint and take the continuum limit. 
For S = 1/2, on the contrary, the scattering length as is 
typically of the order of a few lattice constants and only 
for —4 < J < —3.9 does as become comparable to the 
smallest value of the scattering length for S = 1. 

For S > 1, close to saturation, and in the low-energy 
limit, we therefore map our system onto a dilute two- 
component gas of bosons interacting with effective short- 
range interactions. Many-body effects will be incorpo- 
rated by properly shiftin g th e two-body T matrix off- 
shell as explained in Ref. [32|. We show that, while the 
interspecies interaction is always repulsive and stronger 
than the mtraspecies interaction, the latter undergoes a 
sign change. When the intraspecies interaction becomes 
negative, the bosons are unstable against a collapse. We 
show that a 1/5* expansion captures this physics correctly 
and, similar to the case of the (unfrustrated) Heisenberg 
chain,— is applicable to the present problem, albeit its 
one-dimensional nature. 



vides equivalent results) 

S* = S - ojoj , 5+ = V2Sai , 

S7 = V2Sal(l~al ai /2S), (15) 

we map Eq. (JT|) onto a bosonic problem: 



T (q;k,k') t t 

a k+q a k ,_ q a k a k , , 



k,k f ,q 



21 



where 



(16) 



£ k = Jcosfc + cos2fc — (Jcosfc c ; + cos2fc c ;) > 

is the single- magnon dispersion and k c i = arccos (— J/4). 
Note that in our normalization, the minima of the single- 
particle dispersion are at zero energy: e± kcl = 0. The 
bare interaction vertex Tq is given by To(q;k,k / ) = 
V q - \{V k + V k >) with V k = 2 J cos k + 2cos2fc. The 
chemical potential is /i = hf — h, where hf is the clas- 
sical saturation field value hf = S(J + 4) 2 /4. We are 
interested in the dilute regime (i — > 0. 

Concentrating on the low-energy behavior we arrive at, 
via a Bogoliubov procedure , 25 i 26 a two-component Bose 
gas interacting via a (^-potential with Hamiltonian den- 
sity 

V^ Q | 2 , go(S) 



He 



E 



2m 



-(n 1 +n 2 )-|-go(<S')nin2 . (17) 



Here tp a , a = 1,2 describe bosonic modes with mo- 
menta close to ik c i (i.e., the Fourier transforms of il> a are 
■0i(fc-K)) w a kcl+k , -02(fc->O) w a- kal+k ) and n a = ip^ipa 
are the corresponding densities. The bare coupling con- 
stants of the effective ID model of the two-component 
Bose gas, go(S) and go{S), are, in the dilute limit of 
bosons, related to the renormalized vertices of the micro- 
scopic model Eq. (fTB) through 



r(0; k c i,k c i) = 



T(0; kd,-k c i) 



9o(S) 



r(-2fc ci ; k c i-k c i)-- 



(18) 



go(S) 



i 



go(S)V2m 



The relations Eq. ([18]) follow from a generalization of the 
corresponding equation for the case of a one-component 
Bose ga o 32 ' 36 to the two-component case using an RG 
analysis 3 ^ (see Appendix lAl for details of the calculation). 



B. 1/5 expansion 



A. Effective Hamiltonian 

Using the Dyson-Maleev transformation 3 ^ (the Dyson- 
Maleev representation is used here for convenience. We 
have checked that the explicitly hermitean Holstein- 
Primakoff representation^ to leading order 1/S, pro- 



Next, we apply a 1/S expansion to calculate the inter- 
action vertices and extract the coupling constants go(S) 
and go(S). Using a standard ladder approximation the 
Bcthe-Salpeter equation for the vertices T reads: 

T(q;k,k') = T (q;k,k') 



2SL ^ 

p 



F ( g- p;k+p 7 k' -p) 

^~k-\-p T" £fe' —p 



(19) 
T{p;k,k'). 
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Setting the transferred momentum q = in T and the 
incoming momenta to k = k' = k c i , we get (see Appendix 
E] for details): 



T(0;k d ,k c i) 



T{p;kci,k c i) 



-y 



25L 



(Vb - Vfc cl ) [T(p; fc e ;,fc e ;) -r(0;fc c ;,A: e i)] 

e fcd+P "r e k c i—p 



(20) 



Now, in the spirit of the 1/5 expansion, we replace the 
renormalized vertex with the bare vertex on the right 
hand side of Eq. ([20)) , T(p; k c i,k c i) — > 2e p , which is possi- 
ble since there are no infrared divergences. Regularizing 
the left hand side of Eq. (|2T)]) as in Refs. [3^l3rj and using 
Eq. (HHJ) we extract the coupling constants of the effective 
model. The analytical expression for go(S) is 



9o(S) 



F 



1 



f(J 2 -8) 
|J|S(16-,/ 2 ) 3 / 2 



(21) 



(the derivation of Eq. ([21]) can be found in Appendix [A] 
and the constant F is given in Eq. ([Al2]> ). To leading 
order in J + 4, 



lim g (S) 

,7^-4+ 



S S c 
45 



(J+4) 2 + o((J + 4) 5 / 2 ) 



(22) 



To first order in 1/5, we obtain S cr = 6, which is not 
that large a number, hence corrections beyond 1/5 may 
affect S cr . In the same way we calculate go{S) as 



9o(S) = 



1 



J 2 -8 
16S 



(23) 



(the derivation is presented in detail in Appendix [A] sec 
Eq. (|A25[) for the expression for F) . 

From Eqs. ([21]) and ([23]). we notice that g (S) < g {S) 
for —4 < J < 0. Thus the state below saturation is 
a single-component one. Provided the interactions are 
repulsive, the ground state is a translationally invariant 
chiral stated where bosons prefer to 'condense' at the 
same minimum of the single-particle dispersion since they 
experience a minimal repulsion there We also note that 
\go\ m "C 1 for J < 0, hence interactions between bosons 
are generically weak. In particular, even though m — > 
oo when approaching the ferromagnetic Lifshitz point, 
\g \m 0. 

The effective bare intraspecies interaction is depicted 
in Fig. [5] and behaves as go(S) ~ [ J — J cr (S)} for J — > 
J cr (S). The scattering length is related to the effective 
coupling constant by Eq. ([T4")) . signaling a resonance at 
J cr (5). Thus, we see that for 5 < S cr there is a finite 
region near J ~ —4 where go(S) < and bosons attract 
each other, producing a collapsed state. 



2 10 



IT) 
b/j 




-2 10 



FIG. 2: (Color online) Effective bare intraspecies interaction 
go(S), for S = 3 (solid line), representative of the generic 
behavior for S < S cr , S = 10 (dashed curve), representative 
of S 1 > Scr [dot-dashed line: go(>5 = oo)]. Inset: go for 5 = 1 
and 3/2. 



To corroborate this, using ED for Eq. ([1} and 5=1 
with periodic boundary conditions, we have calculated 
the ground-state momentum of the states with a small, 
but finite number of magnons, which is incommensurate, 
supporting the picture of a uniform chiral state in the 
repulsive case J > J cri and a collapsed state in the at- 
tractive case J < J cr at one of the two minima of the 
single-particle dispersion. 

In the attractive case, d 2 Eo/dn 2 < 0, where Eq and 
n are the bosons' ground-state energy and density, re- 
spectively. In the language of spins, the inverse magnetic 
susceptibility at saturation becomes negative, and hence, 
following standard arguments^ we conclude that there 
is a first-order transition at M = 1, i.e., a jump in the 
magnetization curve just below saturation. 

As pointed out above, the case of 5 = 1/2 is special 
since the scattering length is typically of the order of the 
lattice constant here. The mapping of the 5=1/2 case to 
a two-component Bose gas (by the procedure presented 
above for 5 > 1/2) can be trusted only for J — >• —4, 
where the scattering length becomes much larger than 
the lattice constant. In that case we can easily incorpo- 
rate the exact hard-core constraint into our formalism 2 - 
and again expect that 5 = 1/2 also shows metamag- 
netic behavior. This conclusion is in agreement with 
DMRG results for 5 = 1/2 Note that a metamagnctic 
jump can also be stabilized for spin 1/2 with suitable 
anisotropic exchange interactions However, with our 
procedure we cannot account for the formation of stable 
two-, three-, and four-magnon bound states that is char- 
acteristic for most of the region J > —4 in the spin-1/2 
frustrated ferromagnetic Hcisenberg chain^ 2 - -14 ' 20 

Going back to 5 > 1/2, at lower M. corresponding to 
higher densities of magnons, the hard-core nature of spins 
eventually prevails as well, resulting in a uniform ground 
state at a nonzero momentum. However, as already men- 
tioned, from the finite-size analysis of the two-magnon 
problem, we observe that bound states disappear with de- 
creasing L, suggesting that the attractive effective poten- 
tial (in the limit of a small magnon density) can become 
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(a) S=1,L=128 




0.5 1 1.5 2 2.5 0.07 0.08 0.09 0.1 
h h 

FIG. 3: (Color online) (a) Magnetization curves M(h) for 
S = 1 at J = -2.5, -3, -3.5. (b) M(h) for S = 1, J = -1 
(c) M(h) for S = 3/2 at J = -3.5 (all for L = 128). 



repulsive upon increasing the magnon densities. Thus, 
the state below the jump [i.e., < M < 1 — AMj ump , 
where AMj ump is the height of the jump) will be similar 
to the one encountered in the case of J > J cr , i.e., it is 
a translationally uniform chiral state. 



IV. DMRG RESULTS 

Next, we turn to numerical results for the case of 5 = 1 
(unless stated otherwise) , solving for the ground state of 
Eq. ([T]) in a finite magnetic field h, using ED where pos- 
sible or DMRG.— We present data from DMRG simu- 
lations using up to 1200 states, for L < 128 sites, and 
for open boundary conditions (OBC), unless stated oth- 
erwise. 



A. Magnetization curves and magnetization 
profiles 

The main result of this work, namely, the metamag- 
netic transition from a gapless finite-field phase to full 
saturation, is clearly seen in the magnetization curves 
shown in Fig. [3ja). For 5 = 1, we observe the appear- 
ance of this jump for —4 < J < —2, while for 5 = 3/2 
[an example is shown in Fig. EIc)] the jump exists in a 
much narrower window J < — 3. This is consistent with 
our analytical results for J cr , listed in Tab.|TJ Moreover, 
for — 2 < J < 0, we resolve a plateau in M(h) at M = 0, 
which is due to the Haldane gapi£ in this 5=1 system 
[see Fig. [31(b)] . This gap defines the critical field h c that 



S=1,L=128 



J=-3 (jump at S =97) 

J=-l (no jump) 




FIG. 4: (Color online) Magnetization profiles for S = 1, L = 
128 at J = — 1 (no jump, dashed lines) and J = —3 (jump at 
S z = 97) for S z = 80, 90, 110, 120 (top to bottom). 



separates the gapped Haldane phase from the finite-M 
phase with a smooth M(/i)-behavior for h c < h < h sa f 

The collapse of magnons manifests itself in the magne- 
tization profiles (1 — (5?) vs. site i) using OBC, displayed 
in Fig. [4] for J = —3 and J = — 1. In the former case, 
there is a jump, but in the latter, there is none. Clearly, 
in the states that get skipped over (5 Z > 97 for J = —3), 
magnons collapse into the center of the system, whereas 
for actual ground states below the metamagnetic transi- 
tion, the magnetization profiles become flat. By contrast, 
in the case of J = — 1 where the transition to the fully 
polarized state is smooth and continuous, all profiles are, 
apart from boundary effects, flat. 



B. Central charge 

Our effective theory developed in Section IIHI suggests 
that the gapless phase in the region h c < h < /i sat is a 
one-component phase (where /i sat is the saturation field). 
To substantiate this result, one can make use of entan- 
glement measures such as the von-Neumann entropy to 
extract the central charge, which directly yields the num- 
ber of components of the gapless state. 

The von-Neumann entropy is defined as 



S v n(1) = -tr(/9( In ^) . 



(24) 



where pi is the reduced density matrix of a subsystem of 
length I of our one-dimensional chain of length L. In a 
gapless state that is conformally invariant, the I and L 
dependence of the von-Neumann entropy is given b y 44 ' 45 



S v n{1) = % hi ( — sin(^Z) 
3 \ir L 



9, 



(25) 



which is valid for systems with periodic boundary con- 
ditions (PBC). PBC are preferable for the calculation of 
the central charge from Eq. ([25")) since for OBC, there 
may be additional oscillatory terms, g is a non-universal 
constant that depends on M. As DMRG directly accesses 
the eigenvalues of these reduced density matrices^ it is 
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Spin S=l, PBC 
(a) J=-l,M=l/2 



(b) J=-2, M=l/2 




□ DMRG 
— Fit with c~l 



(c) J=-3, M=5/16 
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FIG. 5: (Color online) DMRG results for the von-Neumann 
entropy S v n{1) in the gapless phase h c < h < h sat of the 5 = 1 
system: (a) J = -1, M = 1/2, (b) J = -2, M = 1/2, (c) 
J = -3, M = 5/16 (symbols). The lines are fits to Eq. ((231) . 
resulting in c = 1.0 ± 0.1 in all cases (we exclude S v n(1) for 
/ < 10 and / > 54 from the fit). In this figure, we display 
results for periodic boundary conditions and L = 64 sites. 
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FIG. 6: (Color online) Phase diagram of Eq. (TT]) for S = 1 
(circles: saturation field /i sa t; squares: spin gap h c ). Inset: 
height AMjump of the metamagnetic jump vs. J. The com- 
parison of L — 64 (open symbols) and L = 128 (solid symbols) 
as well as finite-size scaling (not shown here) supports that 
both AMjump and h c are finite in extended regions of J. 



straightforward to measure S v n(1) with this numerical 
method. 

Some typical DMRG results (squares) for systems of 
L = 64 and periodic boundary conditions are presented 
in Fig. [5J We have fitted the expression Eq. ([2"5j) to our 
numerical data (shown as solid lines in the figure) and ob- 
tain c = 1.0 ± 0.1 in all examples. Therefore, we expect 
the gapless phase to be a (chiral) one-component liquid. 
Note, though, that at both small M and \J\, where the 
convergence of DMRG is notoriously difficult, we can- 
not completely rule out the presence of a c = 2 region, 
which, however, is irrelevant for the main conclusions of 
our work. 



C. Phase diagram for S = 1 

Our results for the 5 = 1 chain are summarized in the 
h vs. J phase diagram Fig. [5j We identify three phases: 
(i) a gapped M = phase at h < h c (similar to the 
Doublc-Haldane phase known for J > 0, see Ref . l4fj|) . (ii) 
a gapless (chiral) finite-field phase for h c < h < h sa ,t , and 
(iii) the fully polarized state at ft, sat < h (with /i sat = 
for J < —4). 

AMj um p is plotted in the inset of Fig.[5j the jump sets 
in at J < —2 (close to where the zero- field gap becomes 
small rendering it difficult to resolve it numerically) , con- 
sistent with our theory. 

Since in the limit of J — > 0, one has two spin-1 
chains with antiferromagnetic interactions which both 
separately have a Haldane gap at zero field ^ upon 
coupling the chains, one obtains the so-called Doublc- 
Haldane phase (in contrast to the regular Haldane phase 
that is inherited from a single spin-1 chain with antifer- 
romagnetic interactions). Both phases, Double-Haldane 
and Haldane phase, are realized in the frustrated, anti- 
ferromagnetic spin-1 chain,— yet in our case, only the 
Doublc-Haldane phase exists. The determination of the 



corresponding spin gap h c is a bit subtle. Namely, a Hal- 
dane chain with open boundaries gives rise to spin-1/2 
excitations at the open ends4£ Since we have two chains 
(for small | J|), we have a total of four spin-1/2 end spins. 
Hence, the spin gap in Fig. [5] is determined from 

h c = E{S Z = 3) - E(S Z = 2) , (26) 

where E{S Z ) is the ground-state energy in a sector with 
a given total S z . 

It is worth emphasizing several differences with the 
phase diagram of the spin-1/2 version of Eq. ([T]). First, 
for S = 1, there are no multipolar phases, which oc- 
cupy a large portion of the corresponding spin-1/2 phase 
diagramA ^ 14 ' 20 Second, the spin-1/2 system features an 
instability towards nematic order ; 11 ' 12 ' 18 which can be 
excluded on general grounds for integer spin- S chains^ 
even for < M < 1. 



V. SUMMARY 

In conclusion, wc showed that resonances can play a 
crucial role in determining the low-energy behavior of 
frustrated quantum spin systems subject to a magnetic 
field. The proximity of resonances caused by an interplay 
between frustration (J > —4) and quantum fluctuations 
(1/2 < S < S cr ) results in extremely large values of the 
ID scattering length that allows to develop an effective 
theory of a weakly interacting two-component Bose gas. 
The quasi-collapse of the dilute gas of magnons provides 
the physical origin of the emergent metamagnctism. 

The preditictions of our analytic theory were verified 
by numerical data. We focused on the case 5 = 1 since 
there the jump in the magnetization curve is the most 
pronounced. As a by-product, we obtained the phase 
diagram for the S=\J-J' chain with antiferromagnetic 
J' and ferromagnetic J in a magnetic field. This phase 
diagram is remarkably simple. In particular, there are no 
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multipolar phases in the case S 
to the spin- 1/2 caseAi^^S 



1, in marked contrast and 
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Appendix A: Details on the mapping to a dilute 
Bose gas 

In Sec. IIII1 we mapped the spin Hamiltonian Eq. ([1]) 
close to saturation to an effective, bosonic field theory 
Eq. (jTTJ) . Here we provide the details of generalizing the 
procedure of obtaining the coupling constants from the 
many-body T-matrix for an effective single-component 
Bose gas model, which is described in Ref. [H, to the 
case relevant to us here, namely an effective theory of a 
two-component Bose gas. 

We introduce the 'mean-field interaction coefficients', 

g(S) = T(0;k ch k cl ) (Al) 
g(S) = r(0;k cl ,-k c i)+T(-2k cV ,k cl ,-k c i), (A2) 

where T(q;k,k') is the full interaction vertex of the 
Hamiltonian Eq. (TIB]) . They satisfy mean- field like re- 
lations for go > go, 



and for g < g 



fi = g(S)n 



V = (g(S)+g(S))n/2, 



(A3) 



(A4) 



where n = (J2 a n a) — ^ is the total density of bosons. 
To get the connection between g(S) and g(S) on the one 
hand and the bare coupling constants of the effective ID 
model of a two-component Bose gas, go{S) and go(S) on 
the other hand, we generalize the corresponding equation 
for the case of a one-component Bose gas^ to the two- 
component case, consistent with the SU (2) symmetry of 
the RG fixed point of a dilute, two-component Bose gas^i 
lim^o g(S) = lim^o g(S) = ny/Ji/y/2m, 



g(S) 



So (SO 



l+.go(S)V2W(7V^)' 
9o(S) 



l + g (S)V2m/(Tr^/fJ.) 
To zeroth order in 1/S 1 , we have 

g (S = oc) = Vo-V kel = ^^ 



(A5) 
(A6) 

(A7) 



g(S = oo) = V 2kcl +V - 2V kal > g{S = oo) . (A8) 

Thus, classically (S —> oo,m — » 0), we have go > go for 
—4 < J < 0. We will show below that incorporating 
quantum fluctuations does not modify this relation. 

In the following we will calculate the effective in- 
traspecies and interspecies interaction constants, to see 
how 1/5 corrections modify ga(S) and go{S)- In order to 
obtain g(S) we have to set in T(q, k, k') given by Eq. (|19[) . 
q = and k = k' = k c i, 



g(S) 



V - V k 
1 



to 



2SL 



2SL 



p 

Vo + Vo 



V kal 



e k c i+p + e fe ci -p 



(A9) 



where we have denoted T(p; k c i, k c i) = T p . After straight- 
forward manipulations we obtain, 



g(S) 



Va - V Kl 



E 



^ £fc c i+P + e k al ~p 



Vo - V kal + 



1 V- 



2SL 



1 - 



V P -Vo 



e k c l+p + £ fcel— P 



-y 



2SL 



(Vo~V k J(T p ~To) 

e fed+P ^k c i—p 



(A10) 



Now, on the right hand side of Eq. (|A10|) . we plug 
in the zeroth order vertex (in the 1/5* expansion), 
r(p; k c i, k c i) 2e p (which is possible because there are 
no infrared divergences anymore) , and use that V p — Vq = 
2(e p — eo). Equation (|A10|1 . after passing to infinite sys- 
tem size, becomes 



g(S) 



i 



i 



4tt5 J tk cl +P + e K,-p 



= F . 



(All) 



where 

F = V -(l 

Vo - V kc 



2S 



2ttS 



)V kcl 



dp 



1 

^5 



dp 



e p( e P ~ e o) 
e kd+p + £fc ci - 



e k c i+p + e fc c (— p 



(A12) 



and we have used that Vo — V kcl = F+0(1/S). According 
to the scheme that we follow the two-body T matrix must 
be calculated off-shell, thus the denominator in Eq. (jAllj) 
must be understood as e p — > e p + Cfi/4:S, where the exact 
value of the numerical constant is C = 7T 2 /8M This leads 
to (compare Eqs. (HHJ) and (|A"T|) ) 



g(S) 



9o(S) 



l + go(S)V2^/(n^)' 



(A13) 
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where we have introduced the intraspecies Licb-Linigcr 
coupling constant go(S) as in Eq. (f2"Tj) . Note that all 
integrals presented in this section are evaluated analyti- 
cally, which is a nice feature of the 1/S treatment of our 
problem. 

Now we outline the calculation of g(S). We denote 
r(p; k c i, — k c i) = T p and obtain the Bethe-Salpeter equa- 
tion for 

TT TT 

f = Vo- V kcl + J^[dpf p - 1 ^[ dp V - p ~ Vkcl f E 



AnS 



and 



h = 



^ ^ l dp V E +V- 2Kl - p -2V kcl 



AttS 



2c 



-kd 

TT — k c 



(f P - f ) 



r f dp Vp + V-™^-^ . (A20) 



AttS J 1 2e kcl+p + Cfi/2S 

-feci 



and 



T-2A; ci = V-2k cl — Vk cl + 



AnS 



dpT f 



f dp V - 2k ; l : p ~ Vk °' t p . (Ai4) 



AttS 



2e k , 



i+p 



Adding these two equations gives g(S) = T + T_ 



2/,-,. 



g(S) = V + V- 2kcl - 2V Kl + 



AttS 



dpl\ 



We divide the last term in Eq. (|A15|) into two pieces, 

7T 

r dp v P + v-2 kcl - p -2v kcl f ^ = h+hj (A16) 



AttS 
where 



2e&, 



i+P 



AirS 



and 



— 7T — k c i 

x(f p -f_ 2kal +t_ 2kel ) 



_^ V f l dp V P + V-2 kel - P -2V kel 



(A17) 



47r5 7 2e fccl+J3 

-feci 

x(f p -f + f ). (A18) 

Note that, for convenience, we have shifted the first Bril- 
louin zone, (— n, tt) — > (— tt — fc c ;,7r — k c i)- Shifting the 
T-matrix off-shell, we have (the integral with a dash de- 
notes the principal value) 

-k c i 



h = — 



r dp ^+ V -^- P -2V kc 



AirS 



— 7T — k c 



2£fe c; 



(f p -f_ afc J 



-p 



AttS J aP 2e kcl+p + C^/2S ^ 

— TT — kd 



The last terms in 1\ can be written as 



7' dp Yp+y-2k cl - P -2v kcl 



AttS J ' 2e kcl+p + Cfi/2S 

-TT-kd 

f-2fe o! dp F p + F - 2fce ,- P -(Fo + F- 2fcgi ) 



4ttS* 



2e 



-7T — fe c 



kd+P 



^-2 kc i f dp v i+ v - 2k r~Z k o- (A2i) 



4^5 



-7T — fc c 



2e fec!+p + C///2S 



Similarly, for the last terms in I 2 we get, 



* f cl dp y P + v-2kci- P -2v kcl 



AttS J r 2e kcl+p + C[i/2S 

-kd 

. . v p + v_ 2kc ^ p -(v + v_ 2kcl ) 



— i 

4ttS J 



dp- 



2e fc , 



■7T—k c 



l+P 



f dp V Q + V-2 kc i - 2V kc 



AttS J ' 2e kel+p + Cfi/2S 



(A22) 



Noting that 



TT — k c 



' dp V P + V -2Kl-p-(V°+V-2kci) 



AnS 



2e kc 



i+p 



-kd 



AttS 



dp 



V p + V- 2kcl -p-(V + V- 2kel ) 



-7T — /c c 



2efe ci 



J — 8 
165 ' 



(A23) 
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and gathering all contributions, Eq. (|A15|) takes the form 

dp 



_ 8 v + V _ 2kci -2V kcl f 
165 16tt5 J 



V a + V 2kcl - 2V kel + 



J 



45 



A^sf p 2Z (T p . kcl -T. 2kcl ) 



r dp Vp - k " + ^+ fc " ^(f p _ fcci -f ).(A24) 



4tt5 



2e„ 



Now we can plug the zeroth-order vertices into the 
right hand side of Eq. (|A24j) . in the spirit of the 1/5 
expansion: T p — > 2e p . Noting that Vo + V-ik al — 2Vfc c! = 
F + 0(1/5), where 



F = Vo + V 2Kl - 2V Kl + 



J 2 



45 



(A25) 



1 



2tt5 
1 

2^5 



d P— f e (ep-fcc -e-2fc OI ) 



An V P -k cl + V P +k cl - 2Vfc ci 

dp Y e \ e v-Ki - eoj, 



in the same way as for g(S) we obtain 



~g(S) 



3o (5) 



l+.9o(5)V2^/(7r v ^) 



(A26) 



with g (S) given in Eq. (|23[). 

We see that in order 1/5, quantum fluctuations do not 
modify the relation go(S) > go(S) for J < 0. Therefore, 
the interspecies interaction is always positive (repulsive) 
for any J and 5 > 1/2, and behaves as <?o(5) ~ ( J + 4) 2 
for J — > —4 for all 5. A hard-core constraint, stemming 
from the mapping of spins to bosons, is not included in 
our theory. While it is well-known how to treat the hard- 
core constraint in an exact way for 5 = 1/2 j2£ this is not 
the case for 5 > 1/2. 



Appendix B: DMRG results for periodic boundary 
conditions 



The DMRG results presented in the main text were 
computed for open boundary conditions (OBC), except 
for those in Fig. [5] where we used periodic boundary con- 
ditions (PBC). Alternatively, one may compute all the 
data, e.g., the magnetization curves M(h) using PBC. 
This approach is generally expected to suffer from (i) 
slower convergence with respect to the number of states 
kept in the DMRG runs2£ and (ii) large finite-size ef- 
fects due to the incommensurability in the problem. We 




-0.5' 



0.5 1 1.5 2 2.5 3 
h 



FIG. 7: (Color online) Magnetization curves M(h) for S = 1 
(a) J = —3.5 and (b) J = — 1 (L = 48), calculated with 
DMRG on systems with periodic boundary conditions. In 
panel (a), we clearly see the metamagnetic jump of height 
AMj ump and in panel (b), the Haldane gap that defines h c 
shows up as a zero-field magnetization plateau. 



here demonstrate that all main features can be seen with 
both OBC and PBC, namely the existence of the meta- 
magnetic jump and the Haldane gap. Moreover, the 
quantitative results are comparable, except for the ex- 
pected finite-size effects due to the incommensurability. 



Magnetization curves from periodic boundary 
conditions for 5 = 1 



First, we discuss some examples of magnetization 
curves obtained for PBC. We kept up to m = 1600 states 
for the PBC DMRG computations presented here. In 
addition, we used exact diagonalization for those sectors 
which arc sufficiently small. Therefore, in particular the 
data close to the saturation field are free of truncation 
errors. 

Figure [7] shows the magnetization curves for J = —3.5 
[panel (a)] and J = -1 [panel (b)] for L = 48 and PBC. 
As in the results for OBC, we observe the presence of the 
metamagnetic jump (here in the case of J = —3.5) and 
the Haldane gap (see the J = — 1 curve), which manifests 
itself as a plateau in the magnetization curve at M = 0. 
The latter defines the critical field h c that separates the 
gapped Double-Haldane phase from the gapless finite-M 
phase (compare Fig. [5]). 

Note that at intermediate M, the PBC data show spu- 
rious small steps with AS Z > 1. We have checked that 
these features disappear as one goes to larger system 
sizes. 



2. Comparison of OBC vs PBC for S = 1 

Figure [5] contains the results for the Haldane gap that 
defines the critical field h c separating the gapped zero- 
field phase from the gapless phase at finite magnetiza- 
tions, the saturation field h sa .t, and the jump height (in- 
set), comparing data for OBC (open symbols) with data 
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J 



FIG. 8: (Color online) Phase diagram of the frustrated ferro- 
magnetic 5=1 chain, comparing DMRG results from systems 
with OBC (open symbols) to results from systems with PBC 
(solid symbols). L = 64 in all cases. Lines are guides to the 
eye. 

from PBC (solid symbols). Here we choose a chain length 



of L = 64 for both PBC and OBC. 

For PBC we can determine the spin gap from 

h% BC = E(S Z = 1) - E(S Z = 0) , (Bl) 

where, as in Eq. (f26|) . E(S Z ) is the ground-state energy in 
a sector with a given total S z . The good agreement be- 
tween the OBC and PBC results for h c in Fig. [8] confirms 
that it is indeed appropriate to use Eq. (|26|) for OBC. 

The PBC data for AM jump in Fig. suffer from the 
presence of several peaks, resulting in a non-monotonic 
dependence on J. This is due to the incommensurabil- 
ity in the finite-magnetization region, which is incom- 
patible with the lattice vectors of a system with peri- 
odic boundary conditions that is translationally invari- 
ant. Therefore, in the main text we focus on the discus- 
sion of DMRG data from systems with OBC. The results 
from OBC and PBC data for the saturation field h sa ,t, 
however, agree very well with each other. 
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